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ABSTRACT 

We present a 0.72 sq. deg. contiguous l.lmm survey in the central area of the COS- 
MOS field carried out to a fcr ~ 1.26 mJy beam"^ depth with the AzTEC camera 
mounted on the 10m Atacama Submillimeter Telescope Experiment (ASTE). We have 
uncovered 189 candidate sources at a signal-to-noise ratio S/N^ 3.5, out of which 129, 
with S/N^ 4, can be considered to have little chance of being spurious (< 2 per cent). 
We present the number counts derived with this survey, which show a significant excess 
of sources when compared to the number counts derived from the ~0.5 sq. deg. area 
sampled at similar depths in the Scuba HAlf Degree Extragalactic Survey (SHADES, 
Austermann et al. 2010). They are, however, consistent with those derived from fields 
that were considered too small to characterize the overall blank-field population. We 
identify differences to be more significant in the Si.imm 5 mJy regime, and demon- 
strate that these excesses in number counts are related to the areas where galaxies 
at redshifts z 1.1 are more densely clustered. The positions of optical-IR galaxies 
in the redshift interval 0.6 ^ z <i 0.75 are the ones that show the strongest correla- 
tion with the positions of the l.lmm bright population (5*1. imm ^ 5 mJy), a result 
which does not depend exclusively on the presence of rich clusters within the survey 
sampled area. The most likely explanation for the observed excess in number counts 
at l.lmm is galaxy-galaxy and galaxy-group lensing at moderate amplification levels, 
that increases in amplitude as one samples larger and larger flux densities. This effect 
should also be detectable in other high redshift populations. 

Key words: surveys - galaxies: evolution - cosmology: miscellaneous - infrared: 
galaxies - submillimetre 



1 INTRODUCTION 

The Cosmological Evolution Survey (COSMOS) 2 sq. deg. 
field has been extensively targeted by a wide array of obser- 
vations in order to probe the cosmic evolution of galaxies and 
the large-scale structure in which they are immersed (Scov- 



ille et al. 2007a). With a wealth of multi- wavelength data 
spanning from X-rays to radio- wavelengths, and a core deep 
UV-optical-IR survey with the highest resolution and sen- 
sitivity offered by space facilities (HST, Spitzer, GALEX), 
it provides a unique opportunity to study the relationships 
and interactions among galaxy populations selected at dif- 
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ferent wavelengths and across a wide array of environments 
in cosmic time. 

A key contribution towards this knowledge comes from 
the far-IR to millimeter wavelength regime, which has 
been shown to uncover ultraluminous violently star-forming 
galaxies at high redshifts [z 2) that would have gone un- 
detected at traditional optical-near IR survey wavelengths 
due to their intrinsic high obscuration (Small, Ivison & Blain 
1997, Barger et al. 1998, Hughes et al. 1998). Named the 
(sub-)millimeter galaxy population (SMG for sort), this pop- 
ulation has been linked to the formation of massive ellip- 
tical galaxies, with large luminosities L ^ IO^'^Lq, large 
star formation rates ^ lOOOM© yr~^, large reservoirs of 
gas ^ 10^°Mq, and large dynamical ^ 1O^^M0 and stel- 
lar masses ^ 10" M0 (e.g. Greve et al. 2005, Tacconi et al. 
2008, Dye et al. 2008). 

Previous mm-wavelength surveys in COSMOS have 
covered areas which were significantly smaller than the full 
2 sq. deg. design survey and the conclusions derived from the 
limited number of detected galaxies suffered from large field- 
to-field variations. COSBO (Bertoldi et al. 2007), a 1.2mm 
survey carried out with the Max-Planck Millimeter Bolome- 
ter Array (MAMBO) mounted on the 30m Institute for Ra- 
dioastronomy at Millimeter Wavelengths (IRAM) telescope, 
detected 37 candidate sources at a S/N^ 3.5 in the central 
0.09 sq. deg. of COSMOS, which was mapped at a Icr noise 
level of ~ 1 mjy beam~^. An adjacent 0.15 sq. deg. field was 
imaged at 1.1mm to a la level of 1.2 to 1.4 mJy beam~^ 
with the AzTEC instrument (Wilson et al. 2008a) mounted 
on the 15m James Clerk Maxwell Telescope (JCMT), uncov- 
ering 50 candidate sources at a S/N^ 3.5 (Scott et al. 2008). 
The combined area sampled by these two surveys amounts 
to only ~ 12 per cent of the total l°.4x 1°.4 COSMOS area. 

Field-to-field variations in the derived overall proper- 
ties of SMGs, such as number counts, redshift distributions, 
or clustering, had already been reported in these kinds of 
smaller (;$ 0.25 sq. deg.) fields and attributed both to vari- 
ance due to the intrinsically small volume sampled by the 
surveys (e.g. Coppin et al. 2006, Aretxaga et al. 2007, Wei/3 
et al. 2009, Austermann et al. 2010) and to the chance ampli- 
fication by a foreground population of galaxies (Almaini et 
al. 2005, Austermann et al. 2009). Indeed, the environments 
traced by optical-IR galax;ies in the COSMOS sub-fields were 
quite different: while the AzTEC survey focused on an area 
with an overdensity of optical-IR galaxies, COSBO sampled 
lower galaxy-density environments (see figure 1). 

A more representative survey of the COSMOS field was 
thus necessary in order to investigate the culprits for these 
large variations and to characterize the global blank-field 
population at 1.1mm. AzTEC alone has surveyed to date 
~ 2 sq. deg. of the blank-field extragalactic sky at 1.1mm, to 
« 0.4 to 1.4 mJy beam~^, coupled to 15 and 10m telescopes 
with resolutions 18 and 30 arcsec, respectively (Scott et al. 
2008, 2010, Perera et al. 2008, Austermann et al. 2009, 2010, 
Hatsukade et al. 2011, Ikarashi et al. 2011). The COSMOS 
survey of 0.72 sq. deg. presented in this paper is the largest 
single-field extragalactic area mapped at 1.1mm at these 
depths, and provides important constraints especially when 
compared to the Scuba HAlf Degree Extragalactic Survey 
(SHADES) fields mapped at 1.1mm, of comparable exten- 
sion and depth (Austermann et al. 2010). 

Section 2 presents the AzTEC/ ASTE observations in 
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Figure 1. Representation of the optically-IR selected galaxy den- 
sity map in the COSMOS field (Scoville et al. 2007b) and the 
uniform coverage areas of the mm surveys published to date: 
the MAMBO COSBO survey at 1.2mm to a Icr noise level of 
~ 1 mJy beam-i (Bertoldi et al. 2007); and the AzTEC/ JCMT 
survey at 1.1mm to a Icr level of 1.2 to 1.4 mJy beam~^ (Scott et 
al. 2008). The AzTEC/ASTE survey at 1.1mm, presented in this 
paper, has an average Icr noise level of 1.26 mJy beam~^ within 
the 0.72 sq. deg. uniform coverage area (solid line). Also repre- 
sented is a concentric area that marks the 25 per cent coverage 
area of the survey, whose noise level increases towards the edges 
of the map. 



the COSMOS field; section 3 details the data reduction pro- 
cesses employed to produce our 1.1mm map; section 4 char- 
acterizes the systematic properties of the map and derives 
the catalog of source candidates from our observations; sec- 
tion 5 compares the AzTEC source catalog to other mm and 
radio-wavelength catalogs in COSMOS; section 6 presents 
the number counts over the 0.72 sq. deg. area, which shows 
an excess over those of the 0.5 sq deg. SHADES field; sec- 
tion 7 explores the origins of the excess in terms of cross- 
correlations with the optical-IR galaxy population; and sec- 
tion 8 discusses the results and summarizes our conclusions. 



2 OBSERVATIONS 

We imaged a 2800 sq. arcmin field centered at right ascension 
RA(J2000.0)=10''00™30.00" and declination Dec(J2000.0) = 
2°14'00.0" with AzTEC mounted on the 10m ASTE (Ezawa 
et al. 2004, 2008), located at 4800m in the Atacama Desert 
of Chile. The survey was carried out from October 20 to 
November 30, 2008, during excellent observing conditions, 
with mean zenith opacity, as reported by the ASTE mon- 
itor, r220GHz — 0.05 and values of r220GHz < 0.06 about 
76 per cent of the time. ASTE was operated remotely with 
the N-COSMOS3 network system (Kamazaki et al. 2005) by 
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observers deployed in San Pedro de Atacama (Chile), Mi- 
taka and Nobeyama (Japan), Amherst (USA) and Tonantz- 
intla (Mexico). A total of 112.6 hrs of on-bolometer time 
was devoted to this field, excluding calibration and point- 
ing observations. The area was sampled in raster mode at 
208 arcsec s~^ along 52 arcmin stripes oriented in azimuth, 
spaced by 1 arcmin steps in elevation. Since the array ori- 
entation is fixed in azimuth and elevation, the scan angle in 
RA-Dec continuously changes due to sky rotation, providing 
excellent cross-linking in the final combined image. A total 
of 203 raster-scan maps were acquired in COSMOS during 
the 2008 observing season, each taking 33 minutes of obser- 
vation in a single passage through the scanning area. Out of 
the 144 nominal bolometers of AzTEC, 117 were operative 
during this season. 

Pointing observations were acquired every half an hour 
to an hour, sandwiching every one to two COSMOS raster- 
maps, depending on observing conditions. The bright QSO 
1055-1-018, Si.imm ~ 2 Jy, was used to measure shifts 
from the standard ASTE pointing model. We acquired 
~ 4 X 4 arcmin^ maps of our pointing target with a con- 
tinuous Lissajous pattern (Scott et al. 2010), and fitted a 
bi-dimensional Gaussian to the resulting map. The calcu- 
lated offset-corrections were not implemented in real time. 
Instead, they were fed into a time series, and a linearly in- 
terpolated offset was applied to the telescope-pointing time 
series of each COSMOS raster-map during the reduction 
process. A total of 270 pointing observations were obtained 
for the COSMOS field during the campaign. 

AzTEC maps are calibrated using planets as primary 
calibrators. Each night Uranus or Neptune was imaged in 
order to derive the flux conversion factor for each detec- 
tor. In a single observation of a fleld the typical statistical 
calibration error is found to be 6 — 13 per cent (Wilson et 
al. 2008a). This estimate is supported by the observations 
of PKS 0537—411, where we report a 6 per cent scatter in 
31 single observations (Wilson et al. 2008b). The COSMOS 
data presented in this paper is the weighted sum of 203 
observations taken over 27 days. If the nightly calibration 
uncertainty is as high as 13 per cent, in the absence of sys- 
tematic errors, combining the 27 days of observations leaves 
us with a statistical calibration uncertainty of 2.5 per cent. 
Adding this in quadrature with the 5 per cent uncertainty 
in the brightness temperature of Uranus at 1.1 mm (Grif- 
fin & Orton 1993) gives a 5.6 per cent overall calibration 
uncertainty. 



observations of the QSO 10554-018. With this correction in 
place, the bolometer signals are flux-calibrated and binned 
into 3x3 arcsec'^ pixels. Performing this process for the 203 
observations of the fleld results in 203 nearly independent 
maps which are then co-added to make a preliminary image 
of the sky. 

As in previous AzTEC analyses, we also produce 100 
noise-only realizations of the COSMOS field by jack-kniflng 
the time-stream data on a row-by-row basis as described in 
Scott et al. (2008). These noise maps are used extensively 
in the characterization of the map properties. 

As a deviation from previous AzTEC analyses, we revise 
the technique presented in Scott et al. (2008), used to esti- 
mate the transfer function of the correlated noise removal al- 
gorithm on point sources in the data. The previous technique 
estimated the transfer function by creating simulated data 
that contain only a point source at the center of the fleld and 
subtracting the eigenmodes identified for removal in the raw 
time-stream data. The revised technique differs in that the 
transfer function is estimated by executing the cleaning al- 
gorithm on data that combines the raw time-stream data 
with several simulated point sources distributed throughout 
the field. In the latter case, the eigenmodes identified for re- 
moval are recalculated from the combination of all sources of 
signal, noise and simulated point sources. In detail, several 
Gaussian point sources are added to the time-stream data, 
all map-making steps save optimal filtering are performed, 
and the resulting maps are differenced from the original true 
maps of the sky. The regions near the simulated sources are 
normalized, stacked and rotationally averaged to produce 
the point source 'kernel.' As demonstrated in Downes et al. 
(2011), this technique provides a better estimation of the 
effects of the non-linear correlated noise removal on a typ- 
ical point source in the map. For the COSMOS field, this 
has the effect of increasing the measured fiux and noise val- 
ues by -1-15 per cent relative to the values obtained from 
the prescription in Scott et al. (2008). We use this kernel to 
optimally filter the maps in an identical manner to that in 
previous AzTEC analyses. 

The final step in the analysis is to use the mean power 
spectral density of the noise maps and the newly estimated 
kernel map to construct an optimal filter for point source de- 
tection. This final set of Wiener filtered maps is composed of 
the filtered signal map, the flltered weight map, and the cor- 
responding signal to noise map. The 100 noise realizations 
are flltered with the same Wiener filter and will collectively 
be referred to as "noise maps" . 



3 DATA REDUCTION 

We reduced the AzTEC data in a manner similar to that 
described in detail by Scott et al. (2008), but with an added 
set of steps to account for nonlinearities in our atmospheric 
cleaning technique. For each of our 203 observations of the 
COSMOS field, the raw time-stream data from the instru- 
ment, which includes both bolometer and pointing data, are 
despiked and then "cleaned" of atmospheric contamination 
in a row-by-row manner using our standard principal com- 
ponent analysis technique (see Scott et al. 2008 for a de- 
scription). An astrometric correction is made to all pointing 
signals in the time stream based on a linear interpolation 
of the pointing offsets measured by the bracketing pointing 



4 I.IMM MAP AND SOURCE CATALOGUE 
4.1 Map 

The inner 0.72 sq. deg. map of COSMOS acquired by 
AzTEC at ASTE is shown in figure (2] This area corresponds 
to a minimum coverage of 50 per cent of the maximum cov- 
erage attained in the map, and has been shown to provide 
excellent noise properties for source extraction and overall 
population analysis (e.g. Scott et al. 2010). The raster scan 
strategy for COSMOS translates into very uniform noise 
properties along this section of the map, ranging from 1.23 
to 1.27 mjy beam~^, while at the extreme edges of the map 
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Figure 2. AzTEC 1.1mm map acquired at ASTE. The circles 
represent the 189 S/N^J 3.5 source candidates that have been 
extracted within the 50 per cent coverage area. For reference, the 
circles have a radius of 30 arcsec. 



(not represented), the noise increases quickly as one moves 
away from the center. The overall noise of the map, char- 
acterized by the combined jack-knifed noise simulations of 
the individual raster-maps, is well represented by a Gaus- 
sian distribution of rms 1.26 mjy beam~^ (see Fig.|3}. This 
value is within the noise range of our previous 0.15 sq. deg. 
COSMOS map (Scott et al. 2008). 



4.2 Astrometry 

The astrometric accuracy of the AzTEC map is verified by 
stacking at the positions of 1471 radio sources in the field 
that are not located within 30 arcseconds of ±3.5(t or greater 
peaks in the AzTEC map. The radio source locations are 
taken from the Very Large Array (VLA) 1.4GHz deep mosaic 
of COSMOS (Schinnerer et al. 2010) which has an rms noise 
level of 15/iJy beam~^ in the central 50' x 50' region, and po- 
sitional accuracy better than 1 arcsec. The stacked AzTEC 
map has a bright, PSF-shaped peak with S/N= 13 and an 
offset from the center of the stacked image of 0.6 ± 1.3 arc- 
sec in RA and —2.6 ± 1.3 arcsec in Dec. Both of these mean 
offset values are small compared to the 3 arcsecond pixel 
size of the map and so we consider them to be too small to 
warrant any correction. 



4 r 



o 2 



. 1 1 1 1 1 1 




r //' 




//' 




/ / 


'\\ 
\\ 


/ /■ 


\ \ 


r / / 
/ / 

p ^ 


'\ \ 


■\ \ 


\ \ 


r f i: 


'■). \ 


f '■' 




f 1- 


'A IIPtUji 




■A nin 






1 : 


\\ 1 




•.\ 


: 


'• \ 

1 ■. , 1 , 


-5 

Sl.lmm/Pixel 


5 1 
(mJy/beam) 



Figure 3. Histogram of flux density values in the AzTEC/ ASTE 
map. The thick solid line represents the values of the signal map 
over the 50 per cent uniform coverage area. Signal from astro- 
nomical sources produces both positive and negative pixels due 
to the fact that the AzTEC map has a mean of zero (Scott et al. 
2008). The dashed line represents the average distribution of flux 
densities found in 100 jack-knifed noise realizations of the field. 
This distribution is well approximated by a Gaussian function of 
mean zero and sigma 1.26 mJy beam~^ (dotted line), which we 
adopt as the average noise level in the uniform coverage area of 
the map. 



source candidates are marked with circles in figure [2] and 
listed in Table [TJ together with their measured S/N, flux 
densities, la photometric errors, and deboosted flux densi- 
ties (see § [6)l . All sources in the COSMOS catalog appear to 
be unresolved. 



4.3 Source Catalogue 

The source extraction algorithm employed to derive the can- 
didate source catalog is identical to that used in Scott et al. 
(2008). We identify point sources in the 1.1mm S/N map by 
searching for local ma:xima within 15 arcsec of pixels with 
S/N^ 3.5 inside the 50 per cent coverage region. The 189 
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Table 1. AzTEC/ASTE source catalog and ancillary sub-mm to radio photometry. The columns give (1) AzTEC identification number, (2) source name, (3) S/N of the detection in 
the AzTEC map, (4) measured 1.1mm flux density, (5) deboosted flux density and 68 per cent confidence interval, (6) flux density at 890^m from SMA observations (Younger et al. 
2007, 2009), (7) flux density at 1.2mm from MAMBO (Bertoldi et al. 2007), (8) flux density at 1.4GHz from VLA (Schinnerer et al. 2010), (9) probability of chance association as per 
P-statistics, (10) angular distance to radio counterpart, and (11) Notes on individual objects, including names for the source in other (sub-)mm catalogs, where AzTEC_J refers to the 
AzTEC/JCMT catalog (Scott et al. 2008), AzTEC denotes interferometric SMA detections at 890/im (Younger et al. 2007, 2009), MMJ denotes detections at 1.2mm by the COSBO 
survey performed with MAMBO (Bertoldi et al. 2007), and B indicates detection at 1.1mm by the Bolocam survey carried out at the Caltech Submillimeter Observatory (CSO; J. 
Aguirre, private communication; for a description of the survey see Aguirre et al. 2006) . Also included in these notes are claims for additional radio associations (RA) and spectroscopic 
redshifts, where references are (a) Younger et al. 2007; (b) Younger et al. 2009; (c) Capak et al. 2011, Riechers et al. 2010; (d) Smolcic et al. 2011; and (e) Capalc et al. in preparation 
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133.0 ± 12.0 


0.001 


1.2 




100025.15+022602.2 


5.4 


6.5 ± 1.2 


5-5±i:l 






56.0 ± 11.0 


0.020 


6.5 


AzTEC.J100025. 23+022608.0, B 


100147.43+022455.9 


5.3 


6.3 ±1.2 


5.4±i:i 














100012.96+020124.1 


5.2 


6.4 ± 1.2 






5.2 ±2.0 


56.0 ± 11.0 


0.020 


6.5 


MMJ100012+020125 


100026.82+023132.4 


5.2 


6.3 ± 1.2 


5.3±J:i 






67.0 ± 12.0 


0.043 


9.5 


AzTEC.J100026.68+023128.1 



AzTEC/Cl 
AzTEC/C2 
AzTEC/C3 

AzTEC/C4 
AzTEC/C5 

AzTEC/C6 

AzTEC/C7 

AzTEC/C8 

AzTEC/C9 

AzTEC/ClO 

AzTEC/Cll 

AzTEC/C12 

AzTEC/C13 

AzTEC/C14 

AzTEC/C15 

AzTEC/C16 

AzTEC/C17 

AzTEC/C18 

AzTEC/C19 

AzTEC/C20 

AzTEC/C21 

AzTEC/C22 

AzTEC/C23 

AzTEC/C24 
AzTEC/C25 
AzTEC/C26 
AzTEC/C27 
AzTEC/C28 
AzTEC/C29 
AzTEC/C30 
AzTEC/C31 
AzTEC/C32 
AzTEC/C33 



Table 1. (continuation) 



AzTEC ID 


Source name 
(MMJ+) 


S/N 


'S'l.lmm 

(m) 

mJy 


'S'l.lmm 

(db) 

mJy 


■SsgO^im 5l.2mm 

mJy nijy 


'S'l.4GHz 


^'l.4GHz 


9 
)) 


Notes 


Az 1 hiO / Co4 


100007.77+021151. D 


5.2 


6.3 ± 1.2 


5.3±i-^ 


5.7 ± 0.9 








ATTVTTinnnn^ i noii a r\ 

MMJ 100007+021149 






c; 1 
O.i 


1 1 
D.z ± i.z 


5.2±j J 




fto n 1 ion 
02. U ± iz.U 


0.U41 


y.o 




AzTEC/C36 


095840.29+020514.7 


5.1 


8.6 + 1.7 


6.8±l:| 




168.0 ± 15.0 


0.000 


0.2 






100121. o2+02ol21.o 


5.1 


6.1 ± 1.2 


c 1-1-1-1 




CO n 1 n 

52.0 ± 0.0 


0.031 


8.0 




AZ 1 £j'^ / '^oo 


1 nnnoQ c;Q_Lnooi/iQ o 
iUUUZo. oo-rUZZi^o.z 


o.i 


O.Z ZIl i . z 


K 1-1-1.2 




/IQ n -1- 1 1 n 
4o.U it ii.U 


U.Uzo 


1 .z 




Az 1 n/L^/L^oy 


1 nni Oft Q/i 1 nonnriQ q 


K 1 
o.i 


O.i ± 1.2 


tr 1 4-1.1 












Az i rjU/L/4U 


nnc;nQ/i Td. i noi nov ft 

uyoyo4. /D+U2iy2 / .d 


5.0 


6.1 + 1.2 


c 1 4-1.1 




ctA n _i_ 1 1 n 
O4.0 ± 11.0 


0.063 


11.6 




AzirL/0/041 


100145. 00+O22 129.0 


5.0 


6.0 ±1.2 












A TTTTi/^ T1 nnm n i nooon^ n A fTTTTi/^F o n^fftl 

AzlEO-J 100019.73+023206.0, AzlEC5, z = 3.97^^ 


Az i / 042 


100019.75+023203.4 


4.9 


5.9 ± 1.2 


4.8±j:; 


9.3 ± 1.1 


1 o/"* n 1 1 r* n 

126.0 ± 15.0 

n 1 1 n 

85.0 ± 15.0 


0.000 
0.028 


0.9 
7.7 


Azl iliC/C4o 


1 nr\nr\o co i nonon^* i 

1U0U03. 00+020200.1 


4.9 


6.0 ± 1.2 


A Oil. 2 

4-8±i.i 




A {\ 1 1 1 n 

84.0 ± 11.0 


0.032 


8.2 




AzTEC/C44 


100033.80+014900.2 


4.9 


6.2 ± 1.3 


r: n 1 i 2 

5.o±;:i 












Az i lijl_'/U40 


10UOUo.o4+02o2o7. 1 


A n 
4.9 


6.0 ±1.2 


A 1 1 - 1 












Az i 11/0/040 


10U114.ol+02ooli.y 


4.9 


6.0 ±1.2 


A o_l1 2 

4.8± J 1 




1 00 n 1 ion 

122.0 ± 12.0 


0.021 


6.7 




Az i rjU/04 ( 


nnc;n/ii 10 1 noninc; ft 

uyoy4i.io+u2uiuo.D 


4.9 


5.9 ± 1.2 






000 n _L 1 1 n 
222. U it 11.0 


0.038 


9.0 




AzirjO/04o 


1UUU39. 24+023847. 9 


4.9 


6.1 ± 1.2 


4.9±}:i 




a n 1 ion 

00.0 ± lo.O 


0.003 


2.5 


A —TTCi/"' T1 nnnoo to i nooo /too 

Az 1 ejO-J 1000o6.72+02o54o.o 


Azi Ji(0/049 


1001 Jl. 83+015403. 3 


4.9 


6.5 ± 1.3 


r: o_i_l 2 

5.3±{:| 












Az i rjLy/OoU 


Uy5yoo.lo+U2Uooo.2 


4.8 


5.9 ± 1.2 


A 1 1 . 1 




c;*? n _J_ inn 
57.0 ± 10.0 


0.001 


1.1 




AzlEC/Col 


^ r\r\r\ A r\ in i mr'noo ^ 

100040.19+015923.7 


4.8 


5.8 ± 1.2 














A fTlT7i/~1 //"irri 

AzlEC/C52 


1 nm r I? f^o 1 nf^oi n/* 

100150.23+022106.3 


4.8 


5.7 ± 1.2 


4.7±J:i 




linn 1 "1 n 

119.0 ± 15.0 


0.026 


7.3 




Azl hiC/Ooo 


100122.65+021211.8 


4.8 


5.7 ± 1.2 


A ft_i_i i 

4.6±;:; 












Az 1 rijO/Oo4 


100125.09+015744.0 


A Q 
4.0 


5.7 ± 1.2 


A ftJ-l.l 












AzlbO/Ooo 


r\r\r\r\K "in i nTr^r'on n 

l()()()()5.iy+()1552().() 


4.8 


5.8 ± 1.2 


4-7±i:2 




F n 1 1 "1 n 

65.0 ± 11.0 


0.074 


12.5 




Az i rL/C/ OoD 


uyoyU4.yo+U221o4.D 


4. / 


5.8 ± 1.2 


A ^ 1 1-1 












A r^TT?/^ I r^^l 

PsJL i n/U/Oo r 


Uyoyoo.o4+021o24.4 


4.7 


5.7 ± 1.2 


A fi4-1.2 




Q n _L ion 
52.0 ± 12.0 


0.049 


10.2 




Az i li/O/Ooo 


innnon at i ni /ic;nn ft 
1UUU2U.4 ( +U14oUU.D 


4.7 


7.3 + 1.6 


c; ftj_i.4 




c;n n J_ ion 
59.0 ± lo.O 


0.019 


6.3 




Az i hjO/Ooy 


"innnon Ar\ 1 nnTV"io 

lUUU3U.4U+U237i2.2 


4.7 


5.6 ± 1.2 


4.6+j 2 




1 £j 1 n 1 1 c c\ 

161.0 ± 15.0 


0.018 


6.1 




AzihjO/ODU 


1 nni no on 1 nooi on t 

100128.29+022129.7 


4.6 


5.5 ±1.2 


4.4±J;l 












AzIxIjO/OdI 


1 nm 1 c\ CO i noo a At~i f\ 

100119. d8+02o442.0 


4.6 


5.7 ± 1.2 


A ft_Li.2 

4.6±i 1 




inorrn 1 1 nnn 

10350 ± 1000 


0.018 


6.1 


extended 


Az 1 iIjU/Odz 


10UlUO.o7+02o7oD.o 


A ft 

4.D 


5.9 ± 1.3 


A '7J_1.2 

4.7±i 2 












J\Z 1 HjU/ODO 


095920.75+020111.0 


A ft 

4.D 


6.1 ± 1.3 


A 1 1 ■ 1 












AzTEC/C64 


100139.64+022345.2 


4.6 


5.5 ± 1.2 














AziEO/Ooo 


nnr"n,io oo \ nooi o/^ i 

095943.22+022136.1 


4.6 


5.6 ± 1.2 


A A^\ 2 

4.4±i:f 




1 r" n t ion 

153.0 ± 12.0 


0.046 


9.8 




AzirjO/Obb 


innincT nn 1 noo/^00 '7 

10UlU5.U0+022o32.7 


4.6 


5.5 ± 1.2 






oc n 1 1 1 n 

86.0 ± 11.0 


0.015 


5.5 




Az i rjO/Ob / 


10Ullo.Dl+02Uy41. ( 


4.6 


5.5 ± 1.2 


A 1 1 ■ 1 












Az 1 rjL^ / L^Do 


1 nni 00 /ID 1 nonvQQ k 
lUUiZz.4o+UZU i 00.0 


A ft 
4.0 


r r 1 -1 rt 
0.0 ± 1.2 


A QJ_1.1 












AzTEC/C69 


100138.07+020908.9 


4.6 


5.5 ± 1.2 


4.3±ii 












AzTEC/C70 


100026.04+020314.9 


4.5 


5.5 ±1.2 


4.3±i;j 




547.0 ± 13.0 
49.0 ± 13.0 


0.022 
0.037 


6.8 
8.8 




AzTEC/C71 


095953.83+021847.9 


4.5 


5.5 ±1.2 


4.3±i:J 


5.7± 1.3 


79.0 ± 11.0 


0.017 


6.0 


MMJ095953+021851, B 


AzTEC/C72 


100159.82+020459.8 


4.5 


5.6 ± 1.2 














AzTEC/C73 


100103.62+022856.8 


4.5 


5.4 ±1.2 


4-2±i:; 












AzTEC/C74 


100105.64+022139.6 


4.5 


5.4 ±1.2 


4-2±l:; 




91.0 ± 12.0 


0.025 


7.3 





Table 1. (continuation) 



AzTEC ID 


Source name 
(MMJ+) 


S/N 


'S'l.lmm 

(m) 


"S'l.lmm 

(db) 


Ssgo^im 


Si. 2mm 


Sl.4GHz 


-Pl.4GHz 


e 


Notes 








mjy 


mJy 


mJy 


mJy 






5) 






lUUUo^.cio+UzUlUo.l 


4.4 


5.3 ± 1.2 


4-2±i:2 














Az 1 rliC'/ / D 


iUUUio.zl+UzlzU / -U 


4.4 


F /I 1 1 

5.4 ± 1.2 


A OJ_l.l 














Azl hiV^/CTT 


095900. UU+UioToD. 7 


A A 

4.4 


C /I _!_ 1 

5.4 ± 1.2 


A OJ_1.2 






ct.(\ n _l_ 1 1 n 
09.0 ± 11.0 


0.014 


5.4 






uyoyuz.7o-i-ui5y4z.u 


/! A 

4.4 


6.5 ± 1.5 


A q_l1-4 














Az i /o / y 


riciKciA o i'r\ 1 no 1 Q /I o o 
Uyoy4o. /U+Uzio4o.o 


4.4 


5.4 ± 1.2 


A 0_L1.1 














Azi hjO/OoU 


T nnno o 0/i i noocrro 


4.4 


5.3 ± 1.2 


A 1 _1_1.1 














AzirjO/L8i 


iUUUUo. 97+1)15241. 4 


4.4 


5.4 ± 1.2 


A 1 _l1.2 














Azi rjO/Ooz 


T nn OCT 1 noi(^/io /i 


4.4 


5.2 ± 1.2 


A 1 _l1 














Azi hjO/Ooo 


T nnoon ocr i noi a 1 a f 

lUUzoU.zo+Uzi4i4.D 


4.3 


6.8 ± 1.6 


4.9±J:| 














Azl hjC/Co4 


095942. 7D+U15511. 7 


4.3 


5.3 ± 1.2 


A ^ Jul 1 














Azl hiC/Coo 


1 nr\i .;<n nn i nofvc d i o 

100140. 09+022541.0 


4.3 


5.2 ± 1.2 


4.0+;:} 






Ac\ r\ 1 ion 

549.0 ± 12.0 


0.039 


9.0 




Az 1 HiV^/UoD 


100109.00+021720.0 


4.0 


5.2 i 1.2 


A n 1 1 . 1 






/loin _i_ Af\r\ 
4210 lb 400 


0.016 


5.7 


extended 


Az 1 HjU/Uo f 


1 nnonc i noi '7nn i 
100205.50+021700. 1 


4.0 


5.2 i 1.2 


A n 1 1 . 1 

4.0+1 2 














Azi hjC/Ooo 


nntrnoT ot p non/ioo n 

Uy5yo7.o7+U2U42ci.y 


4.3 


5.3 ± 1.2 


4.0+1 1 














Az i rjO/Coy 


T nni ov nc i noi Qo^i i 
1UU12 i .Uo+U21oo0.1 


4.3 


5.1 + 1.2 


A n_Ll. 1 
4.0+1 1 














Az i rLO / CyU 


inniQc^ i-i A 1 noii^c;n o 
lUUlo0.u4+U21oo0.2 


4.3 


5.1 + 1.2 


A nj-l.l 
4.0+1 2 














Azi hjC/091 


100128.61+022347.4 


4.2 


5.1 ± 1.2 


3.8+1;? 






81.0 + 12.0 


0.038 


8.9 




Azi hjO/Oyz 


100139.94+023015.0 


4.2 


5.2 ± 1.2 


4-o±l:2 














Azi hjC/Cyo 


1 nni m no i noi i OT n 

1001o2. 02+021 lo7.0 


4.2 


5.1 ± 1.2 


3.8+J:; 






i2f\ n 1 1 1 n 

60.0 + 11.0 


0.004 


2.8 




Az 1 to / oy4 


nncnc7 i q j noi '7i n n 
095957. 1O+021719.0 


A 

4.2 


C 1 _!_ 1 

5.1 ± 1.2 


3.9+1 2 














Az 1 CiU/v-^yo 


1 nnm Q oti i noio/io n 
lOOOlo.OD+021242.9 


A 

4.2 


5.1 + 1.2 


3.8+1 1 














Az 1 HiU/ uyb 


1 nm no qo i nicic/t a 
lOOlOo. 00+015154. D 


A 

4.2 


5.2 i 1.2 


3.9+1 2 














Az i HiV^/uy / 


1 nnoi A Q A 1 r\m a A A n 
100214.04+02 1944. y 


A 

4.2 


5.5 + 1.3 


A 1-1-1.2 
o_l1. 1 

3.8+1 2 






55.0 + 16.0 


0.013 


5.1 




Az 1 to/oyo 


1 r\r\r\A o nn i nonc;iQ /i 
100042. yy+02051o. 4 


/I 
4.2 


5.0 + 1.2 






TQ n -L A A n 
78.0 + 14.0 


0.005 


3.1 




Az i / oyy 


t r\r\r\r\(i nc i ni c^nc^o n 

looooo.yo+uioyoo.y 


4.2 


5.1 + 1.2 


oil. 1 

3.8+^ 2 














Azi hjO/OlUU 


nncrnio tro i noinocr o 

095918. 52+021035. S 


4.2 


5.1 + 1.2 


Q 01 1.2 

3.8+1 1 














Azi hjO/OlUl 


rtnc n A r:: cro i noonio r" 

095945.58+020018.5 


4.2 


5.1 + 1.2 


3.8+1:1 














Az i HjO/ ClUz 


0y5y48.77+02ol5o.5 


4.1 


5.1 + 1.2 


1 1 . 1 

3.8+1 2 














Azi tC/ClOo 


in/^iOj< 1 mpT/^i^ 1 

100124.43+015615.1 


4.1 


5.0 ± 1.2 


3.8+1:1 






51.0 + 15.0 


0.036 


8.7 




AzTEC/C104 


095940.74+015333.1 


4.1 


5.3 ± 1.3 


A r\ \ ^ 1 

4.o+{:^ 














AzlEC/ClOo 


095845.11+021442.1 


4.1 


5.5 ± 1.3 


/I 1 _i_l 2 














AzIEC/CIUd 


100006.54+023838.6 


4.1 


5.3 ± 1.3 


4-o±l:i 


8.6 + 1.3 










Az 1 EC_J 100006.40+023839.8, Az 1 EC6 


AziE0/01U7 


r\(\TZ(\r>(\ Jz c 1 nooooo 

0959o9.5o+0222o8.o 


4.1 


5.0 ± 1.2 


3.8+i'2 














Az i rjU / UlUo 


T nm T nrc i nooi^i A q 
1UU11d.Uo+U2od14.o 


A 1 

4. 1 


5.3 + 1.3 


A n-l-1.2 














AzTEC/C109 


100111.63+022838.3 


4.1 


5.0 + 1.2 


3.7+11 






59.0+ 11.0 


0.004 


2.7 




AzTEC/CllO 


100108.56+020029.7 


4.1 


5.0 + 1.2 


3.7+1:1 














AzTEC/Clll 


095929.62+021241.6 


4.1 


5.0 + 1.2 


3.7+i:i 






67.0+ 12.0 


0.019 


6.3 




AzTEC/C112 


100010.94+015309.3 


4.1 


5.0 + 1.2 


3.7+1:1 






122.0 + 12.0 


0.014 


5.3 




AzTEC/C113 


095914.96+022957.7 


4.1 


5.5 + 1.3 


4.0+1:1 






173.0 + 16.0 


0.036 


8.7 




AzTEC/C114 


100024.06+022000.6 


4.1 


4.9 + 1.2 


3.7+1:1 














AzTEC/C115 


100014.84+020532.5 


4.1 


4.9 + 1.2 


3.7+1:1 















Table 1. (continuation) 



AzTEC ID 


Source name 
(MMJ+) 


S/N 


'S'l.lmm 

(m) 

mjy 


•S'l.lmm 

(db) 

mJy 


SsgOfim Sl.2min 

mJy nijy 


'S'l.4GHz 


.Pl.4GHz 




" 


Notes 


Azl hjO/Cllb 


100109.63+020348.3 


4.0 


4.9 ± 1.2 


3.7±i:^ 


4.2 ± 1.4 


59.0 ± 11.0 


0.008 


4.0 


MMJ 100109+020346 


Az 1 r!iC'/v^il / 


oyoy2o.yo+022oio.o 


4.0 


5.0 lb 1.2 


Q '7_l_ 1 . 1 












AzTEC/Cll8 


095959.49+020633.1 


4.0 


4.9 ± 1.2 


3.7+1:1 


5.6 + 1.1 


104.0 ± 13.0 


0.023 


6.9 


MMJlOOOOO+020634, B 


AzTEC/C119 


095915. 34+020748. 3 


4.0 


5.0 di 1.2 


3.7+1:1 












Az i rjO/OlzU 


lUUiU0.oo+U2U2i4.D 


4.0 


4.9 + 1.2 


3.6+1 1 












Az i hjC/Olzl 


095952.00+020915.9 


4.0 


4.9 ± 1.2 


3.7+1:1 












Azi bC/C122 


100200.88+021648.0 


4.0 


4.8 ± 1.2 














AzirL/0/012o 


1UUU22.d2+U1o145.4 


4.0 


4.9 ± 1.2 


3.7+1:1 












AziiljC/L'lz4 


095940. o2H-U2c>o54.0 


4.0 


5.1 lb 1.3 


Q '7_l_1.2 

3.7+1 2 




orr n 1 ion 

27.0 ± 13.0 

1 oncrn 1 1 onn 

loU59 ± loUU 


0.038 
0.033 


8.9 
8.3 


A rrn 771 /~\ '\ r\f\T1i C\ A C Cd 1 nOQC/11 n 

Az 1 l!jO_JU9594D.Db+U2o541.9 
extended 


Az i H^U/Uizo 


U9o92u.o4-f-U22DOo.o 


A n 
4.0 


4.9 lb 1.2 


3.6+1 2 












Az i CjKj/KjLZo 


iUUio9.4o-+-U2zzoy.4 


A n 
4.0 


4.8 lb 1.2 


q c_i_1.2 
O.Oitj 2 




T A n _L 1 Q n 
lo4.U ± lo.U 


0.048 


10.1 




Az i hiLj/LjlZ i 


innioc Act 1 noQr^T/i q 
iUU12o.4oH-U2oo24.o 


A n 
4.0 


5.4 zb 1.3 


Q q_l1!3 
o.o±J 




1 Q 1 n _L ion 
lol.U ± 12. U 


0.007 


3.9 




Az i rLO/dzo 


1 nnnno on i ni c;i 100 
lUUUUz.oU+Ulolio.o 


4.0 


4.9 + 1.2 


3.6+1 2 












A r,nPT7^^ on 
Az i rL/C/dzy 


1 nni on ou 1 nonoi 7 


4.0 


4.7 + 1.2 


Q cr 1 1.1 

3.5+1 1 










ts 


Azi rjO/OloU 


nncrncrr' nn 1 nonono /I 

095956.99+020308.4 


3.9 


4.8 ± 1.2 


3.5+1 2 




1 /I n 1 inn 

124. U ± lU.U 


0.018 


6.1 




Azi rjC/Olol 


"innnin no 1 nooTcr-i 

lUUUlU.9o+Uzo754.d 


3.9 


4.9 ± 1.2 


3.5+1 2 












Az i rjC/doZ 


1 nnnncc qct i noovc;'/ /I 
lUUUUo.oo+U2o ( / .4 


3.9 


4.9 + 1.3 


3.5 + 1 2 












Az 1 rLiLy / ioo 


1 nni on ft/i 1 nooftO/i a 


0.9 


4. / it L.Z 


3.4+1 2 




'7C1 n 1 ion 
(9.U it iz.U 


U.UUO 


Q A 
0.4 




Az ± ll/'^ / io4 


1 nni n^i ni _i_ni r;ni a 
iUUiUO.Ui-rUioUi4.z 


o.y 
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Figure 4. Expected fraction of false detections in the AzTEC 
catalog as a function of limiting S /N estimated through counting 
peaks in 100 noise-only realizations of the field. 



4.4 False detection rate 

The nature of false-positives in source identification in 
AzTEC maps is richly discussed in Perera et al. (2008) and 
Scott et al. (2010). We follow the same procedure as Scott 
et al. (2010) and search our 100 noise realizations of the 
COSMOS field in order to derive a conservative upper limit 
to the fraction of sources in Table [1] that are actually noise 
peaks rather than real sources. Figure U shows this upper 
limit of the fraction of false detections as a function of S/N. 
The catalog is robust. At S/N^ 3.5 we derive a mean of 17.4 
false detections (9 per cent of the source candidates), and 
at S/NJS 4.0 we expect only 3 (2 per cent) to be false. 



4.5 Completeness 

We estimate the completeness of source detections through 
a set of simulations in which artificial sources of different 
flux densities are inserted one at a time within the 50 per 
cent coverage of the observed COSMOS map, and then re- 
trieved with the same source extraction algorithm as that 
used for building the catalog. As described in detail by Scott 
et al. (2008, 2010) this method has the advantage of taking 
into account the effects of random and confusion noise in the 
signal map while it does not alter the properties of the real 
map significantly. The simulations are based on 1000 test 
sources per represented flux value. A test source is consid- 
ered recovered if it is extracted with S/N^ 3.5 at a radius 
^ 17 arcsec of its input position. This radius is adopted to 
ensure a ~ 100 per cent recovery of S/N= 3.5 sources (see 
§4.6). Figure [S] summarizes our findings. Sources with flux 
densities SLimm 5S 6 mjy are identifled in this survey with 
^ 90 per cent completeness. 



Figure 5. Survey completeness estimates for the COSMOS cat- 
alog as a function of flux density. The points and 68 per cent 
binomial error bars show the completeness estimated by inserting 
sources of known flux density one at a time into the observed sig- 
nal map and then finding them with the same source extraction 
algorithm as that used to create the candidate source catalog. In 
order for a source to be considered recovered, it must be detected 
with S/N^ 3.5 at a distance r ^ 17" from the input location. 

4.6 Positional uncertainty 

The large beam of ASTE combined with low S/N detec- 
tions and confusion noise contribute to introducing random 
offsets in the position in which a source is found in a map. 
We characterize the positional uncertainty as a function of 
S /N from the same set of simulations as those used for the 
completeness calculation, where now we focus on the dis- 
tance at which the sources get recovered. Figure [S] shows 
the uncertainty distribution for sources found in three S/N 
regimes, and the comparison with the analytical approxima- 
tion of uncertainties (Ivison et al. 2007), assuming a FWHM 
of 34 arcsec for the effective beam of the map; the prob- 
ability of flnding a source at distance > D from its true 
position is given by P(> D) = 1 — exp(— D^/2a"'^) , where 
(J = Aq = A5 0.6 FWHM (S/N)-^ The analytical form 
shows the same trend as the empirical distribution, and can 
be used to estimate search radii of possible counterparts of 
the SMGs presented in this paper. 

5 COMPARISON WITH OTHER MAPS 

5.1 Overlap with other mm surveys 

The AzTEC/ ASTE map includes areas of the sky previously 
imaged by other mm-wavelength surveys, including COSBO 
(Bertoldi et al. 2007) and AzTEC/ JCMT (Scott et al. 2008). 
Figure[7]represents the catalogs of these two surveys overlaid 
on the AzTEC/ASTE map, in order to emphasize common 
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Figure 6. Positional uncertainty distributions for AzTEC/ASTE 
source candidates. The data-points and 68 per cent confidence 
intervals show the probability P{> D) that AzTEC sources of 
different S/N will be found outside a radial distance D, as deter- 
mined from simulations. The curves show the expected probabil- 
ities from the simple approximation that takes into a account the 
beam size and the S/N of the detection (Ivison et al. 2007). 



sources and differences among catalogs. Also shown are the 
unpublished sources detected by Bolocam at the CSO (J. 
Aguirre, private communication) that are common to the 
AzTEC catalog. 

The AzTEC/ASTE and AzTEC/ JCMT surveys have 
similar noise properties with rms ~ 1.25 mjy beam~^. The 
10m ASTE image, however, has a resolution of 34 arcsec, 
while the JCMT image has an effective resolution of 18 arc- 
sec, after Wiener filtering. Due to the considerable level of 
incompleteness at low-S/N in both catalogs, not all sources 
are expected to be found in both. Out of the 50 candi- 
date sources extracted from the AzTEC/ JCMT survey at 
S/NJS 3.5, 48 fall within an overlapping region with the 
AzTEC/ASTE survey. In order to find which entries in the 
catalogs match each other, we adopt the positional uncer- 
tainty relation that depends on the effective beam size and 
S/N of the detections (§ 14. 6p . This relationship has been 
shown to work well for both the AzTEC/JCMT (Scott et 
al. 2008) and ASTE data (Fig. [6]). Since a real source will 
suffer from positional uncertainties in both catalogs, we will 
use for each ASTE source a search radius at the 95 per cent 
confidence level of containing the location of the real source, 
and will add in quadrature the 95 per cent confidence po- 
sitional uncertainty radius of a potentially matching JCMT 
source. If the distance between the JCMT and ASTE catalog 
positions are smaller than the resulting search radius, both 
catalog entries will be considered to correspond to the same 
source. Search radii derived in this manner range from 6.4 
to 13.6 arcsec. A candidate source that is found in two inde- 
pendent datasets also increases its individual reliability over 



Figure 8. Deboosted ASTE versus deboosted JCMT fluxes cor- 
responding to the common sources between this paper and the 
revised catalog published in Downes et al. (2011). Error bars rep- 
resent 68 per cent confidence intervals in deboosted fluxes. 

the average false-detection rate corresponding to its nominal 
S/N. Out of the 48 AzTEC/JCMT source candidates that 
fall within the ASTE surveyed area, 16 are common to the 
ASTE catalog, and are listed in Table [T] 

The coincidence between the bright AzTEC/ASTE and 
AzTEC/JCMT sources is remarkable. From the list of the 
15 brightest AzTEC/JCMT sources that were followed up 
and confirmed with SMA (Younger et al. 2007, 2009), only 
AzTEC 13 and 14 are not found in the AzTEC/ASTE cat- 
alog. Both sources have deboosted flux densities Si.imm ~ 
4.4 ± 1.3 mJy (Scott et al. 2008). At these flux densities, 
the completeness estimated for the AzTEC/ASTE map in- 
dicates a ~ 57 per cent chance of recovery. The other two 
sources in this flux range, AzTEC/JCMT 11 and 12, are in- 
deed detected. The level of recovery at lower flux densities 
is lower, as expected by the completeness function. 

Figure [8] represents the deboosted fluxes of the 
AzTEC/ASTE map versus the deboosted fluxes of the 
AzTEC/JCMT revised catalog (Downes et al. 2011). The 
filtering and deboosting methods employed were exactly the 
same for both datasets. The mean ratio between fluxes is 
1.08, with an r.m.s. of 0.44. Folded into these dicrepancies 
are both the increased noise in the new flltering technique, 
blending of fainter sources within the larger ASTE beam 
and the impact of these effects on the deboosting process. 

At 1.2mm the COSBO survey derived a catalog of 27 
candidate sources detected at Si. 2mm 2.2 mJy, out of 
which 12 have AzTEC/ASTE 1.1mm counterparts, follow- 
ing a positional uncertainty analysis similar to that per- 
formed above. The FWHM for the MAMBO map has been 
adopted to be 11 arcsec, and the resulting search radii at 
95 per cent confldence level range between 6.2 and 14.0 arc- 
sec. Common sources are listed in Table [T] The mean de- 
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Figure 7. AzTEC map, and comparison with other mm catalogs. Circles and numbers denote AzTEC sources from this work. Blue 
and red squares denote COSBO (Bertoldi et al. 2007) and AzTEC/JCMT (Scott et al. 2008) sources, with accompanying labels C# and 
A# that refer to their catalog numbers, respectively. The boundaries of the surveys are marked with dashed lines. The AzTEC sources 
from this work also detected by the Bolocam survey are marked with crosses (Aguirre et al., private communication). AzTEC sources 
from this work robustly associated with 1.4GHz sources in this or other works are marked with inner yellow squares. 



© RAS, MNRAS 000, 



AzTEC millimeter survey of the COSMOS field - III. 13 



boosted flux density ratio for these common sources is 
Si. 1mm/ 5'i,2mm — 0.93 lb 0.43. Among the 15 brightest 
COSBO sources, with S'i.2mm ^ 4.4 mjy, COSBO ID num- 
bers 2, 6, 11, 12 and 14 are not formally detected in our 
AzTEC/ ASTE catalog, ahhough they show 2.4 ^ S/N 
3.4 peaks in the AzTEC S/N map. For instance, COSBO 
source 2, Si. 2mm 5.9 mJy (Bertoldi et al. 2007), has 
a probability of detection in the AzTEC/ ASTE map of 
~ 90 per cent (assuming a 0.93 mean flux ratio), and, indeed, 
has a S/N« 3.0 in the AzTEC/ASTE map at the COSBO 
position. 

From the catalog of 19 candidate sources extracted 
by the Bolocam team, 10 are also detected by the deeper 
AzTEC/ASTE survey, and are identified in Table [T] The 
adopted effective FWHM for the Bolocam map in the po- 
sitional uncertainty calculation to match AzTEC/ASTE 
sources is 31 arcsec (Laurent et al. 2005), and the result- 
ing search radii range between 10.9 and 18.7 arcsec. 

5.2 Overlap with the deep COSMOS radio-survey 

The VLA 1.4GHz deep mosaic of COSMOS (Schinnerer et 
al. 2010) also provides an excellent catalog to look for coun- 
terparts of the SMGs presented in this paper. The tight cor- 
relation between radio continuum emission, which is domi- 
nated by synchrotron radiation from supernova remnants, 
and FIR emission dominated by thermal radiation from 
warm dust heated by young stars, in galaxies (Helou, Soifer 
& Rowan-Robinson 1985, Condon 1992, Yun, Reddy & Con- 
don 2001) translates into a large percentage of cross iden- 
tifications among catalogs derived in both frequencies (e.g. 
Ivison et al. 2002, 2007). 

The matching process was carried out within 17 arc- 
sec radii circles of the AzTEC positions, looking for possi- 
ble associations. To quantify the significance of the possible 
associations, we have used the P-statistic (Downes et al. 
1986) which calculates the probability that a radio source 
of a given observed flux density could lie at the observed 
distance from the AzTEC source by chance. Only possible 
radio counterparts with values of P ^ 0.05 are considered 
robust, and are listed in Table[T] together with their distance 
offsets and P- values. 

Out of the 189 AzTEC source candidates, 77 (40 per 
cent) have robustly associated radio counterparts within the 
COSMOS radio catalog, and out of these 77, 7 (10 per 
cent) have a double robust association. These percentages 
are lower than those of other SMC studies ( ~ 60 — 70 per 
cent, e.g. Ivison et al. 2007, Biggs et al. 2011) due to the large 
incompleteness at the 1.4GIIz catalog limit (a cut thresh- 
old is used) and shallower nature of the 1.4GIIz COSMOS 
data. This is also the reason why some previously claimed 
associations (see Table 1) are not identified in this study, as 
they were based on lower threshold detections. 

5.3 Multiwavelength photometry 

Table [1] summarizes the photometry derived from the com- 
parison with other sub-mm to radio-wavelength catalogs. 
The 1.4GIIz fiux densities of the robustly associated coun- 
terparts are one to two order of magnitude lower than 
the deboosted 1.1mm fiux densities, with the exception 
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Figure 9. Millimeter to radio flux density ratio as a function of 
rcdshift for 20 local galaxies used as templates in the derivations 
of photometric redshifts of SMGs (Aretxaga et al. 2007). The grey 
horizontal lines represent the colours of the AzTEC galaxies with 
robustly associated radio counterparts, as in Table 1. The arrow 
at the top left corner represents the typical Scr lower limit for 
SMGs in the COSMOS map with undetected radio counterparts. 



of a few marked extended sources (see below), suggesting 
the mm fiux has a thermal origin rather than the hyper- 
luminous synchrotron dominated blank-field sources iden- 
tified by the South Pole Telescope (Vieira et al 2010) or 
the variable Si. 2mm > 10 mJy flat-spectrum quasars dis- 
covered in MAMBO flelds (Voss et al. 2006). None of the 
bright Si.imm(db) ^ 8 mJy AzTEC sources are associated 
with luminous X-ray objects either (Johnson et al. in prep) , 
suggesting that, at best, they might harvour weak active 
galactic nuclei or be Compton-thick. 

In thermal dominated sources, the millimeter to radio 
flux density ratio can be exploited as a redshift indicator. It 
increases monotonically with redshift, with some degeneracy 
due to the variety of radio synchrotron-slopes and mm dust- 
emissivity indices present in the interstellar medium of those 
local galaxies used to define the relationship (Carilli & Yun 
1999, 2000). Additionally there exists a level of degeneracy 
between the temperature of the dust generating the rest- 
frame FIR luminosity (and hence mm fiux) and the redshift. 
Regardless, by adopting a library of local galaxy templates, 
and accepting the intrinsic dispersion in their SEDs, the 
1.4GHz to 1.1mm fiux-density ratio still provides a crude 
but useful estimation of the redshift (Carilli & Yun 2000, 
Rengarajan & Takeuchi 2001, Aretxaga et al. 2007). This 
indicator becomes relatively insensitive to redshifts beyond 
2 ~ 3, as the 1.1mm filter starts to sample the fiattening 
of the SED towards the rest-frame FIR peak, whilst still 
providing a powerful discriminant between low-redshift (z < 
2) and high-redshift (z > 2) objects. 

Figure [9] represents the colour-redshift diagram of 20 



© RAS, MNRAS 000, 



14 Aretxaga I. et al. 



local galaxies used as templates in the derivations of photo- 
metric redshifts (Aretxaga et al. 2007) and the colours mea- 
sured for the AzTEC sources with radio associations. Except 
for a few exceptions (like AzTEC/C61, a likely radio-loud 
AGN, or CI 24 which is associated with a likely foreground 
radio galaxy at z ~ 0.3) that appear at the bottom of the di- 
agram, the majority of the SMGs have mm-to-radio colours 
indicative of z > 1 systems, in accordance to the expecta- 
tions derived from other (sub-)mm-wavelength redshift dis- 
tributions for the population (Chapman et al. 2003, 2005, 
Aretxaga et al. 2003, 2007, Pope et al. 2005, Valiante et 
al. 2007). A detail analysis of radio and optical-IR counter- 
part associations to these AzTEC sources and their inferred 
redshifts is deferred to an upcoming publication. 



6 NUMBER COUNTS 

We derive estimates for the number density of SMGs as a 
function of flux density, the so-called "number counts" using 
the Bayesian technique originally outlined in Coppin et al. 
(2005, 2006) and used extensively in previous AzTEC pub- 
lications (for example, see Austermann et al. 2009). While 
other techniques commonly used for the extraction of source 
counts from submm/mm wavelength surveys, in particular 
the "P{D)" approach (e.g. Patanchon et al. 2009, Glenn 
et al. 2010), can in principle estimate the counts at fainter 
flux densities (i.e. below the detection limit of the survey), 
these methods are highly dependent on the assumed model, 
and the formally derived errors do not always represent the 
true uncertainty in the counts at faint flux densities (e.g. see 
discussions in Scott et al. 2010 and Glenn et al. 2010). On 
the other hand, with the Bayesian technique, the estimated 
counts are only weakly dependent on the assumed model of 
the prior distribution (Austermann et al. 2009, 2010), and 
the derived error bars more accurately describe the uncer- 
tainty in this estimate. For this reason, we use the Bayesian 
method, and derive the source counts only down to a flux 
density limit of 5(1. 1mm) = 1.5 mjy. 

We parameterize the number counts using a Schechter 
function: 

with N being the number of sources per square degree, S 
being the source flux density at 1.1mm wavelength, and a 
being the power law slope of the faint-end counts. In this 
formalism N^rnjy has a natural interpretation as the number 
of sources per square degree with a flux of 3 mJy. 

The measured flux densities of sources blindly detected 
in the AzTEC map must be corrected for biases resulting 
from the interaction between the Gaussian noise distribu- 
tion in the map and the underlying flux density distribution 
of sources on the sky. We perform this correction by con- 
structing the full posterior flux density distribution for each 
source, taking as a prior the parameters A'smjy = 160 sq. 
deg.~^, S" — 1.3 mJy and a — —2.0, which are consistent 
with those measured in the SHADES fields (Austermann et 
al. 2010). We nevertheless iterate on the adopted prior values 
to guarantee weak dependence on the starting values for the 
end results. For the ^ 3.5a peaks in the map, the full pos- 
terior probabilities are parameterized by their maxima and 



68 per cent confidence levels, and are listed as deboosted 
fiuxes in Table 1. 

Once the full posterior fiux distribution is derived for 
all ^ 2.5(7 peaks in the map, we cut all source candidates 
whose posterior fiux distribution indicates a 5 per cent or 
greater probability of having a negative intrinsic fiux. This 
is likely a very strict cut on our catalog, however, the larger 
beam-size and greater confusion in this survey warrants a 
conservative first approach (Austermann et al. 2009). 

We bin our resulting catalog in 1 mJy fiux density bins, 
correct each bin for the corresponding completeness, and 
calculate dN/dS uncertainties by bootstrap sampling the 
dN/dS probability distribution in each bin 20,000 times, 
taking also into account the error in completeness. Table 2 
lists the resulting bin centers, differential number counts, 
and 68 per cent confidence interval uncertainties and fig- 
ure [10] shows the differential and integrated number counts 
derived for the fuU 0.72 sq. deg. COSMOS field in this man- 
ner. Also plotted are the number counts for all other pub- 
lished AzTEC extragalactic blank-fields (Scott et al. 2008, 
Perera et al. 2008, Austermann et al. 2010, Hatsukade et al. 
2011) following a re-analysis of each of these maps using the 
same technique outlined in section 3, in order to ensure data 
processing effects take no role in the differences found. 

Overall, the counts from different fields show some strik- 
ing variance, especially if one focuses on the two largest, 
COSMOS and SHADES: the 0.72 sq. deg. COSMOS field 
with systematically more sources at all fiux levels than 
the 0.5 sq. deg. (to approximately the same noise level) 
SHADES fields. The result however, is not as extreme as 
that shown in the comparison with the smaller 0.15 sq. deg. 
AzTEC/JCMT COSMOS field. In sections 7 and 8 we wiU 
explore further these differences. 

We derive best-fit parameters for the Schechter func- 
tion that describes the COSMOS/ASTE differential number 
counts by fitting Eq. 1 with a Levenberg-Marquardt least- 
squares algorithm that uses the full data-data covariance 
matrix in the calculation. Our data does not meaning- 
fully constrain the faint end of the counts and thus we fix 
the faint-end power-law index, a, to a value of —2. This will 
also allow a direct comparison to similar fits in the litera- 
ture. The best fit values of iVamjy and S" for COSMOS are 
given in Table 3 and Figure ITT] represents the error contours 
for these parameters. 



7 AZTEC SOURCES VS. LARGE SCALE 
STRUCTURE IN THE COSMOS FIELD 

With the detection of ultra-bright SMGs by AzTEC (Wilson 
et al. 2008b, Ikarashi et al. 2011), MAMBO (Lestrade et al. 
2010), the South Pole Telescope (Vieira et al. 2010) and the 
Herschel satellite (Negrello et al. 2010), there has been new 
attention given to the role that lensing plays in our view 
of the SMC population. Our AzTEC survey of the COS- 
MOS field offers a prime opportunity to inspect the rela- 
tionship between the background SMGs and the foreground 
large scale structure mapped out at optical-IR wavelengths 
over a large redshift span. Indeed, with the 0.15 sq. deg. 
AzTEC/JCMT map of COSMOS Austermann et al. (2009) 
already found a signiflcant spatial correlation between the 
projected foreground [z <^ 1) galaxy population and the 50 
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Figure 10. (Left:) Differential number counts for the COSMOS AzTEC/ASTE field (red) along with the re-analyzed counts from 
previously published AzTEC studies. Because of its large field, the COSMOS AzTEC/ASTE counts are more sensitive than previous 
studies in the 3-15 mJy range. Overplotted (lines) are a number of 1.1mm number count predictions from an array of semi-analytic 
models for galaxy formation, where z for the Rowan-Robinson (2009) models denotes the free parameter that describes the onset of the 
IR-luminous phase. (Right;) Corresponding integrated number counts at 1.1mm. 



Table 2. COSMOS/ASTE difi'erential and integral number 
counts, calculated as described in the text. The difi'erential num- 
ber counts flux bins are 1-mJy wide with effective bin centers 
(first column) weighted according to the assumed prior. 
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Table 3. Best-fit Schechter function parameters to the COSMOS 
AzTEC/ASTE differential number counts with alpha fixed to —2. 
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Figure 11. Confidence limits for N-j^jy and S' when holding 
a = —2. Contours represent the 68 per cent and 95 per cent 
confidence limits. 
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SMG candidates. In this section we will take advantage of a 
much wider range of foreground large scale structures cov- 
ered within the 0.72 sq. deg. AzTEC/ASTE map of COS- 
MOS in order to test if the foreground structure significantly 
impacts our view of the blank-field SMG population and 
whether this might be the likely origin of the difference in 
number counts between COSMOS and SHADES. 

Figure [12] shows the projected density of optical-IR 
galaxies in the AzTEC field with photometric redshift 
«phot ^ 1-1 along with the location of the 129 AzTEC 
sources with S/N ^ 4.0, ^ 2 per cent of which we expect 
to be spurious detections. While there is no apparent cor- 
relation between the positions of both populations, we can 
quantify this impression by comparing the distribution of 
projected densities of optical-IR galaxies within 30 arcsec 
radii circles centered on the AzTEC positions with those of 
30 arcsec radii circles centered on random locations in the 
map (as in Austermann et al. 2009). A Kolmogorov-Smirnov 
(KS) D-test can then be used to ask with what probability 
we can rule out the null hypothesis that the two samples are 
drawn from the same parent distribution. As expected, the 
result, Pks = 68 per cent, is consistent with no strong cor- 
relation between the optical-IR population at z 1.1 and 
the S/N > 4 AzTEC catalog. 

Performing the same test with a flux-cut on the AzTEC 
catalog results in a more significant correlation detection. 
Considering only the 41 S/N ^ 5 sources we find Pks = 
1.8 per cent, which is tentative evidence (at « 2.4cr level) 
that the null hypothesis of no correlation is rejected. Indeed, 
for the 20 AzTEC sources with deboosted fiux-densities 
•S'l.imm ^ 6 mjy we can strongly reject the null hypothe- 
sis of identity between the distributions of galaxy-densities 
around random positions in the AzTEC coverage area of the 
COSMOS field and galaxy-densities around bright AzTEC 
sources: Pks = 0.11 per cent (« 3.3o-). 

The latest release of photometric redshifts for optical- 
IR selected galaxies in the COSMOS field has achieved ac- 
curacies A2/(l + z) ^ 0.007 to 0.012 at 2 ^ 1.25 (Ilbert et 
al. 2009). This statistical precision allows for the identifica- 
tion of large-scale structure pertaining to different redshift 
slices (see for example Scoville et al. 2007b). With this in- 
formation in hand, we can search in redshift space for the 
structures that are more likely associated with the bright 
AzTEC sources. We will again compare the distributions of 
galaxy-densities around AzTEC sources and around random 
positions in the AzTEC covered COSMOS area at the red- 
shift of interest. As a reference, if we take into account the 
AzTEC catalog of 129 S/NJ^ 4 sources, there is no significant 
signal of statistical difi'erences between the distributions at 
any redshift (see Fig. I13p . 

However, if we restrict the analysis to the 20 sources 
with deboosted fiux densities Si.imm 6 mJy, we identify 
the redshift ranges 0.58 ^ z <, 0.76 and 0.21 ^ z ^ 0.26 
as the ones in which most significant (~ 3cr) differences 
are found between the median density of galaxies around 
AzTEC sources and the median density of galaxies around 
random positions (Fig. [13]). The same redshift bins are high- 
lighted if the analysis focuses on the 42 sources with de- 
boosted fiux densities S'l.imm ^ 5 mJy, although the proba- 
bility of rejecting the null hypothesis for the former redshift 
interval is smaller than for the later. 

The optical-IR galaxy density in these redshifts slices 
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Figure 12. Smoothed surface-density map of galaxies with pho- 
tometric redshifts 2 < 1.1 derived from the optical-lR catalog of 
COSMOS (Scoville et al. 2007b), updated to include photometric 
redshifts derived from 30 optical-IR bands (Ilbert et al. 2009). 
Only the 0.72 sq. deg area surveyed by AzTEC with uniform 
(^ 50 per cent) coverage is represented. Darker colours indicate 
more densely populated areas of the sky. The cross and plus sym- 
bols represent AzTEC sources detected at signal-to-noise ratios 
S/N^ 4 and 4 > S/N > 3.5, respectively. The contours divide the 
map into zones with lower or higher than the average density of 
optical-IR galaxies. The total areas of these zones are 0.392 sq. 
deg. and 0.329 sq. deg, respectively. 
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Figure 13. Probability that the median galaxy density around 
AzTEC sources is significantly larger than the density of galaxies 
around random positions in the AzTEC-covered map versus red- 
shift. The thick black solid line represents probabilities around 
deboosted Si.imm 5? 6 mJy sources, the grey solid line represents 
probabilities around deboosted 5i.imm ^ 5 mJy sources, and 
the dashed brown line represents probabilities around S/N ^ 4 
sources. The redshift bins have an increasing step ranging from 
Az = 0.014 to 0.026, at z = 0.15 to 1.1, to accommodate a 
sufficient sample of galaxies and the increasing degradation in 
photometric redshift precission (e.g. Scoville et al. 2007b) 
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Figure 14. Smoothed surface-density map of galaxies at 0.58 <J 
z <j 0.76 detected at optical-IR wavelengths by the COSMOS 
survey within the uniform coverage area of the AzTEC/ ASTE 
map. The cross symbols represent 20 AzTEC sources detected 
with deboosted flux densities SLimm ^ 6 mjy. The plus sym- 
bols denote sources detected with deboosted flux densities 5 ^5 
•Sl.imm < 6 mJy. A massive cluster at z 0.73 (Guzzo et al. 
2007) located in the NW of the map is marked as #1. The inner 
yellow circle of 1.5 arcmin diameter marks the core of X-ray emis- 
sion, while the purple 6 arcmin diameter circle marks the FWHM 
of the optical/IR overdensity associated with the cluster (Scoville 
et al. 2007b). Another rich cluster at z ^ 0.61 from the large-scale 
structure catalog of COSMOS (Scoville et al. 2007b) is marked as 
#24, and the FWHM of the optical/IR overdensity is encircled 
in purple. The outer blue dashed-line contour depicts the edge of 
the uniform coverage of the AzTEC map observed in 2005 at the 
JCMT (Scott et al. 2008). 



and the positioning of the bright AzTEC sources can be 
seen in figures. [14] and [151 In the interval 0.58 ^ z <y 0.76 
we find that Si.imm ^ 6 mJy AzTEC sources have a signif- 
icantly denser galaxy environment than that found at ran- 
dom positions in the map: the null hypothesis of identity 
between the medians of the distributions can be rejected 
at a 1 — Pmw = 99.88 per cent confidence level, while for 
Si.imm ^ 5 mJy AzTEC sources that level gets reduced to 
1 — Pmw = 98.7 per cent. In the interval 0.21 ^ z <^ 0.26 
the reverse happens, finding that the Si.imm 6 mJy 
AzTEC sources have a denser environment with a signifi- 
cance 1 — Pmw = 99.2 per cent. The significance gets in- 
creased to 1 — Pmw = 99.95 per cent if Si.imm ^ 5 mJy 
sources are considered. The effect seems to be carried by 
four 5 ^ S'l.imm < 6 mJy sources that coincide with large 
density peaks, rather than by the general population of in- 
termediate/brightness sources, though. If one excludes the 
S'l.imm ^ 6 mJy sources from the analysis, the probability 
drops to 1 — Pmw ~ 99.0 per cent. 

These two redshift intervals are the same ones we iden- 
tified as having the largest correlation between z <, 1.1 
optical-IR galaxies and AzTEC sources for the smaller 



Figure 15. Smoothed surface-density map of galaxies at 0.21 < 
z <J 0.26 detected at optical-IR wavelengths by the COSMOS sur- 
vey within the uniform coverage area of the AzTEC map observed 
in 2008. Symbols and lines are as in Fig. 1141 Two X-ray clusters 
from the catalog of large-scale structures in COSMOS (Scoville 
et al. 2007b) are marked as #17 and #22. 



0.15 sq. deg. field observed by AzTEC in the JCMT (Auster- 
mann et al. 2009) . It is thus important to know if the corre- 
lations arise mainly due to the very rich clusters located in 
the area of the sky previously sampled, or if this is a trend 
observed over the larger 0.72 sq. deg. field, which might be 
more representative of a generic blank-field. 

In Austermann et al. (2009) we showed that the associ- 
ation of bright AzTEC sources around areas of large galaxy 
densities was not exclusively due to the presence of mas- 
sive rich clusters. If we mask out the areas marked as #1 
and #24 in figure [T!] which represent the FWHM of the 
two massive clusters at 0.58 <, z <, 0.76 identified in the 
large-scale-structure analysis of Scoville et al. (2007b), and 
exclude the sources that fall within those areas, the proba- 
bility of having a median of galaxy-densities around AzTEC 
sources larger than that of random positions by chance is 
only Pmw ~ 0.3 per cent. The distributions are shown in 
figure [TS] Likewise, if we exclude from the correlation anal- 
ysis the overlapping area common to the AzTEC/ JCMT 
and AzTEC/ ASTE maps (encircled in dashed blue lines 
in figures [2] and [TS} , the remaining 12 Si.imm ^ 6 mJy 
sources still show a larger tendency of falling within the large 
galaxy-density regions mapped at optical-IR wavelengths at 
z ^ 1.1. In that case, the null hypothesis of identity be- 
tween the medians of the distributions can be rejected at a 
1 — Pmw = 97.0 per cent confidence level. 
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Figure 16. Histograms of the number of optical/IR galaxies 
around random positions in the AzTEC covered COSMOS map, 
normahzed to the mean value over the whole map, represented a 
solid black line. The dashed, dash-dot-dash and dash-dotted lines 
represent the distribution of number of galaxies found around 
AzTEC sources with deboosted flux densities Si.imm ^ 6 mjy 
in the entire AzTEC covered COSMOS map, excluding the 6 ar- 
cmin diameter area around clusters #1 and #24, and excluding 
the whole uniform area mapped by AzTEC in 2005 in the JCMT. 
The vertical bars at the top of the diagram mark the mean values 
of the distributions. 
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Figure 18. Distribution of fluxes for the 60 brightest (Si.imm ^ 
5 mJy) AzTEC sources that fall within high galaxy-density ar- 
eas (continuous line) and within low galaxy-density areas (dotted 
line) of the COSMOS map. The high density areas host the ma- 
jority of the brightest sources. The null hypothesis that the two 
distributions can be derived from the same parent distribution 
can be rejected with a probability Pks = 0.37 per cent. 



8 DISCUSSION 

8.1 Effects of Foreground Structure on Number 
Counts 

With evidence that our detected SMGs are spatially cor- 
related with foreground {z ^ 1.1) large scale structure in 
the COSMOS field, we now ask what effect, if any, this 
structure might have on our estimation of the SMG number 
counts. We address this by first dividing the catalog into 
two sub-samples: SMGs that lie in regions of the map with 
foreground galaxy densities larger than the average density 
of the field and SMGs that lie in regions with foreground 
galajcy densities smaller than the average galaxy density of 
the field. Figure [12] shows the two regions along with the 
respective locations of the corresponding SMGs. 

Figure [17] shows the number counts resulting from split- 
ting sources in high and low foreground galaxy density envi- 
ronments. While none of the dN/dS bins individually show a 
clear, statistically significant deviation from the mean counts 
over the field, the data clearly suggests that, consistent with 
the results of Section 7, at high fiuxes (Si.imm 6mJy) 
the counts are systematically higher for the high foreground 
galajcy density sample. This result can also be made evident 
by representing the distribution of flux densities of the SMGs 
that fall within the low and high galaxy-density areas of the 
map (see Fig. llSp . The null hypothesis that the 1.1mm flux- 
density distribution of sources that fall within low and high 
density areas are similar can be rejected, as differences as 
large as the one measured can only be produced by chance 
in 0.37 per cent of situations. 

The number counts for the low density galaxy sample 
is closer to the number counts derived for SHADES, if still 



systematically above them (by ~ 16 per cent in flux den- 
sity), while the high density galaxy sample aligns with the 
results derived for COSMOS/ JCMT and GOODS-N, (about 
~ 30 — 40 per cent offset in flux density at Si.imm ^ 5mJy, 
compare with Fig. 10). It is thus apparent that an accurate 
description of the overall population will fall somewhere in 
between these solutions, and that despite having sampled 
0.72 sq. deg. in COSMOS and 0.5 sq. deg in SHADES at 
comparable noise levels, this is not yet a large enough area 
to avoid variance due to intrinsic clustering of the SMG pop- 
ulation and ampliflcation by foreground structures. 

The optical-IR galaxies that compose the COSMOS 
galaxy density map have accurate photometric redshifts that 
place them at z ^ 1.1, while most of the redshifts of the 
SMGs uncovered by AzTEC are still unknown. Furthermore, 
the largest amplitude correlation between the positions of 
optical-IR galaxies and SMGs occurs at 0.6 ^ z <y 0.75 
(§ 7), where only a ~ 3 per cent of 850/im SMGs with ra- 
dio associations are statistically located (Chapman et al. 
2005). Photometric arguments (Fig. 9, § 5.3), place the ma- 
jority of AzTEC SMGs at z ^ 1 (see also Younger et al. 
2007, 2009) and the 4 bright targets with reflned interfero- 
metric positions by SMA that have been spectroscopically 
targeted and have robust redshifts are indeed at z > 1.1 
(Table 1, Capak et al. 2011 and in prep., Riechers et al. 
2010, Smolcic et al. 2011). While there is no proof that all 
bright Si.imm ^ 6 mJy sources are at high redshifts, it is 
likely that most of these sources are at z > 1 and that their 
association with the optical-IR galajcy large scale structure 
is through ampliflcation, an explanation already posed to 
account for the discrepancies between the number counts 
measured in the smaller COSMOS AzTEC/ JCMT field and 
SHADES (Austermann et al. 2009). Our new results, thus. 
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Figure 17. (Left:) Differential number counts derived for regions of low and high galaxy density within the AzTEC/ASTE covered 
COSMOS area. "High foreground" denotes sources that fall in regions of the map with greater than the mean foreground galaxy density. 
"Low foreground" denotes sources that fall in regions of the map with less than the mean foreground galaxy density. Triangles denote 
counts from the COSMOS/JCMT. (Right:) Corresponding integral number counts. 



confirm this interpretation for a more representative area 
of the sky that is not dominated by the presence of a rich 
cluster environment. 

Lensing of galaxies by foreground galaxies and by fore- 
ground groups of galaxies has been shown to be the reason 
for the presence of a very bright (sub-)mm selected galaxy 
population (Vieira et al. 2010, Negrello et al. 2010, Lestrade 
et al. 2010, Ikarashi et al. 2011). These galaxies are extreme 
cases of the phenomena presented here, where moderate 
~ 20 per cent amplification factors would be enough to ac- 
count for the fiux density offset between the number counts 
of SHADES and COSMOS. 

Light propagation experiments in cosmological simu- 
lations using multiple lens-plane algorithms show that at 
z > 1 any population is subject to a large chance of am- 
plification by foreground structures (e.g. Martel & Premadi 

2008) . Correlations between bright SMGs and foreground 
optically selected galaxies at 2: ~ 0.5 were first detected 
using a small sample of S'sso^m ^ 8 mJy sources in the 
UK 8-mJy and Hubble Deep Field (Almaini et al. 2005). 
The analysis of a statistically robust sample of 2477 350/im- 
selected SMGs in the Lockman-SWIRE field, that includes 
the smaller UK 8-mJy survey area studied previously, has 
confirmed correlations between the SMC population and 
z ^ 0.2 and z ~ 0.4 optical and IR-selected galaxy samples 
(Wang et al. 2011). Our first 0.15 sq. deg. survey in the COS- 
MOS field scanned all foreground structures at < 2: < 1, 
yielding significant correlations between the bright 1.1mm- 
selected SMGs detected in the field and optical-IR galax- 
ies at redshifts z ~ 0.25 and z ~ 0.65 (Austermann et al. 

2009) . Whether these correlations were dominated by the 
optical-IR overdensity of galaxies where the AzTEC survey 
had been acquired or inherent of amplifications to the gen- 



eral blank-field population was open to debate. Our new 
analysis in COSMOS allows for a better estimation of the 
structures that contribute to the possible boosting of flux 
densities of SMGs by foreground structures, identifying the 
0.21 <, z<, 0.26 and, most signiflcantly, the 0.58 <, z <, 0.76 
redshift bins as those with the largest probability of associ- 
ation with bright Si.imm ^ 5 mJy sources. The first redshift 
interval is common to the 350/im-selected SMG correlations 
result. The second redshift interval, however, is still unex- 
plored by other experiments. 

Other populations of high-z galaxies, also should show 
similar amplification trends. The positions of bright Lyman 
Break Galaxies selected to be at z ~ 2.5 — 5 in the Canada- 
France-Hawaii- Telescope Legacy Survey are, indeed, posi- 
tively correlated with optically-selected galaxies at z < 1.4, 
and this effect has been shown to be consistent with weak 
lensing by the foreground structures in the line of sight 
(Hildebrandt, Waerbeke & Erben 2009). Strong lensing is 
also predicted to be a dominating effect in the bright num- 
ber counts of SMGs ( Blain 1996, Negrello et al. 2007, Lima, 
Jain & Devlin 2010), however, neither strong nor weak lens- 
ing are often included in the predictions of observables of- 
fered by galaxy formation models. 



8.2 Comparison with galaxy formation models 

Figure 10 shows the number counts of COSMOS compared 
with several semi-analytical galaxy formation models that 
have successfully reproduced overall properties of the SMG 
and SMG population (Granato et al. 2004, Baugh et al. 
2005, Rowan-Robinson 2009). The 1.1mm number counts 
of the SHADES field were compared with these very same 
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models (Austermann et al. 2010), and it was found that they 
all overpredicted the number counts in the > 3 mjy regime. 
The discrepancies with the new 0.72 sq. deg. AzTEC/ASTE 
COSMOS field are not that severe. While the Granato 
et al. (2004) model for the joint formation of QSOs and 
SMGs could be made compatible with the COSMOS num- 
ber counts using a small shift in flux density that would 
mimic the amplification claimed in this paper, the model 
over-predicts the 5i.imm ^ 2 mJy number counts, which, 
although they are poorly constrained by the COSMOS data 
alone, have more robust estimations from smaller deeper 
fields (Scott et al. 2010, Hatsukade et al. 2011). The Baugh 
et al. (2005) model offers a good description of the COS- 
MOS number counts at Si.imm ^ 4 mJy, but it over-predicts 
the number counts at lower flux densities, while the Rowan- 
Robinson (2009) models have a different functional form to 
that found for COSMOS over the whole flux density range 
sampled by our study. 



9 CONCLUSIONS 

The number counts of the COSMOS 0.72 sq. deg. 
AzTEC/ASTE field show an overdensity of sources with re- 
spect to the population of SMGs found in previous large 
blank-field surveys such as SHADES (Austermann et al, 
2010). The number counts derived for the COSMOS field 
display a systematic positive offset over those of SHADES, 
but are nevertheless consistent with those derived from 
smaller fields that were considered too small to characterize 
the overall blank-field population. We identify departures 
to occur more significantly in the 5*1. imm ^ 5 mJy regime, 
and demonstrate that these differences are related to the 
areas where galaxies at redshifts z <, 1.1 are more densely 
clustered. The positions of optical-IR galaxies in the red- 
shift interval 0.60 ^ z ^ 0.75 are the most strongly cor- 
related with the positions of the 1.1mm bright population 
(SLimm ^ 6 mJy), a result which does not depend exclu- 
sively on the presence of rich clusters within the survey sam- 
pled area. The most likely cause of these departures in num- 
ber counts at 1.1mm is lensing by either foreground galaxies 
or foreground groups of galaxies at moderate amplification 
levels, that increases in amplitude as one samples larger and 
larger flux densities. Our results and the comparison with 
the previously published SHADES number counts illustrate 
the fact that even ~ 0.70 sq. deg. surveys are still subject to 
variance due to the small volume sampled by the mapped ar- 
eas in conjunction to the chance amplification by foreground 
structures. 
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